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Relativistic Hartree-Bogoliubov (RHB) theory in coordi- 
nate space is used to describe the chain of even-even Zir- 
conium isotopes reaching from ^^®Zr to the drip line nucleus 
^''"Zr. Pairing correlations are taken into account by a density 
dependent force of zero range. For neutron numbers larger 
than the magic number A'^ = 82 a giant neutron halo out- 
side the core of ^'^^Zr is observed. It is formed by up to six 
neutrons. 

PACS numbers : 21.10.Pc, 21.10.Gw, 21.60.-n,21.60.Jz 



With perpetual improvement in the techniques of pro- 
ducing radioactive beams, the study of exotic nuclei very 
far away from the line of /3-stability has now become fea- 
sible. Experiments of this kind may cast new light on 
nuclear structure and novel and entirely unexpected fea- 
tures may appear: Neutron rich nuclei can have a struc- 
ture very different from that of normal nuclei. They con- 
sist of a normal core surrounded by a skin of neutron 
matter. Close to the drip line, where the coupling to the 
continuum becomes important, a neutron halo can de- 
velop, as it has been observed in several light nuclei, the 
most famous being '^-'^Li 1^. However, in all of the ha- 
los observed so far, one has only a very small number of 
neutrons, namely one or two outside of the normal core. 

In order to study the influence of correlations and 
many-body effects it would be very interesting to find 
also nuclei with a larger number of neutron distributed 
in the halo. In this letter we report on the theoretical 
prediction of a giant neutron halo for Zr-isotopes close to 
the neutron drip line. It is formed by up to six neutrons 
outside of the ^^^Zr core with the magic neutron number 
N = 82. 

Recently, a fully self-consistent relativistic method has 
been developed for the description of halo nuclei. It is 
based on Relativistic Hartree-Bogoliubov (RHB) theory 
in coordinate space It is able to take into account 
at the same time the proper isospin dependence of the 
spin-orbit term, which is very crucial for a reliable de- 
scription nuclei far away from the line of /3-stability ^] 
and the self-consistent treatment of pairing correlations 
in the presence of the continuum As it has been 
shown in Ref. the scattering of Cooper pairs into the 
continuum containing low-lying resonances of small angu- 
lar momentum plays an important role for the formation 
a neutron halo. By using a density dependent zero range 
interaction, the halo in ^-'^Li has been successfully repro- 
duced in this self-consistent picture. Very good agree- 



ment with recent experimental data is obtained without 
any adjustment of parameters. To obtain these results, a 
full solution of the RHB equations Q in coordinate space 
is necessary. The expansion in a harmonic oscillator basis 
[Q , which is very useful only for nuclei close to the line of 
/3-stability, provides only a very poor approximation in 
the continuum, even if a large number of oscillator shells 
is taken into account. The simple BCS-approximation 
used in Ref. [Q leads to a partial occupation of unbound 
states in the continuum and therefore to a gas of evap- 
orating neutrons. Here we report on calculations, where 
the same method as in Ref. has been used to investi- 
gate halos in Zr-nuclei at the neutron drip line. 

The starting point is the Relativistic Mean Field La- 
grangian 



C = ip\^^- g^(/j - QppT - -e(l - T3)4 - g„(j -mjip 
+ ^d^<jd>'a - U{a) - ^n,,M^'' + ^mlijpij^ (1) 

It describes the nucleons with the mass m as Dirac 
spinors tp moving in the fields of mesons: a scalar me- 
son (cr), an isoscalar vector meson (w), an isovector 
vector meson (p) and the photon A^, with the masses 
TO(j, and irip and the coupling constants go-, g^j, 
gp. The field tensors for the vector mesons are given 



as rj^jy — dp 



duUjp and by similar expressions for 



the p-meson and the photon. For simplicity we ne- 
glect these fields in the following, however, they are fully 
taken into account in the calculations. For a realistic de- 
scription of nuclear properties a nonlinear self-coupling 
U{(t) — \m?^a'^ -\- \g2(y^ + iflso''* for the scalar mesons 



has turned out to be crucial 



Using Green's function techniques it has been shown in 
Ref. 1^ how to derive the relativistic Hartree-Bogoliubov 
(RHB) equations from such a Lagrangian: 
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Ek are quasiparticle energies and the coefficients Uk{r) 
and Vfe(r) are four-dimensional Dirac spinors. h is the 
usual Dirac hamiltonian 



(3) 



containing the chemical potential A adjusted to the 
proper particle number and the meson fields cr and uj de- 
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termined as usual in a self-consistent way from the Klein 
Gordon equations in no-sea-approximation. 
The pairing potential A in Eq. (||) is given by 

Aa. = ^EKTc.'^c. (4) 

cd 

It is obtained from the pairing tensor k = U*V^ and the 
one-meson exchange interaction Vjbcd ™- pp-channel. 
More details are given in Ref. |q|. There it has been 
found, that these forces are not able to reproduce even 
in a semi-quantitative way proper pairing in the realistic 
nuclear many-body problem. As in Ref. we therefore 
replace V^^^^ in Eq. (|^) by the density dependent two- 
body force of zero range: 

4 \ po J 

(5) 

where r = \ri + r2|/2. 

The RHB equations (|) for zero range pairing forces 
are a set of four coupled differential equations for the HB 
Dirac spinors U{r) and V{r). They are solved for the 
parameter set NLSH Q , which is widely used for the de- 
scription of neutron rich nuclei, by the shooting method 
in a self-consistent way. The details will be published 
elsewhere. With a step size of 0.1 fm and using proper 
boundary conditions the above equations are solved in a 
spherical box of radius i? > 15 fm. The results do not 
depend on the box size for i? > 20 fm. For a radius 
i? = 30 fm we found the same results within an accu- 
racy of 0.1 % for the binding energies and matter radius. 
Since we use a pairing force of zero range we have to 
limit the number of continuum levels by a cut-off energy. 
For each spin-parity channel 20 radial wavefunctions are 
taken into account, which corresponds for i? = 20 fm 
roughly to a cut-off energy of 120 MeV. For fixed cut-off 
energy and for fixed box radius R the strength Vq of the 
pairing force (|^) for the neutrons is determined by a cal- 
culation in the nucleus ^^^Zr adjusting the corresponding 
pairing energy —^TiAk to that of a RHB-calculation us- 
ing the finite range part of the Gogny force DIS |10[. For 
Pq we use the nuclear matter density 0.152 fm~^Tln or- 
der not to miss any bound state the cut-off energy has 
to be larger than the depth of the potential. But as long 
as this is the case and as long as the interaction strength 
is properly renormalized, the results of this investigation 
stay practically unchanged. And in such a way, we get 
almost the same result from the the pairing force of zero 
range and the finite range Gogny force. 

In the upper panel of Fig. 1 we show the rms radii of 
the protons and neutrons for the Zirconium isotopes with 
mass numbers A = 110 to ^ = 140, the drip line nucleus. 
They are given as 



where A^^ = Z, N for protons and neutrons and 

{rh - (/ r^pi^d^r^ (7) 

are the contribution of the different blocks with the quan- 
tum numbers I and j. To guide the eye we also give a 
dashed line with a A^^/'^-dependence. It clearly shows a 
kink for the neutron rms-radius at the magic neutron 
number N = 82. 

These results can be understood more clearly by con- 
sidering the microscopic structure of the underlying wave 
functions and the single particle energies in the canoni- 
cal basis ||ll[. Therefore, in the lower panel of Fig. 1 we 
show the single particle levels in the canonical basis for 
the isotopes with an even neutron number as a function 
of the mass number. Going from A^ = 70 to A^ = 100 we 
observe a big gap above the I/111/2 orbit. For neutron 
numbers larger than the magic number A^ = 82, the neu- 
trons are filled to the levels in the continuum or weakly 
bound states in the order of 8^3/2, 2/7/2, 3pi/2, 2/5/2 
and lft.9/2- 

The neutron chemical potential is given in this Figure 
by a dashed line. It approaches rapidly the continuum 
already shortly after the magic neutron number A^ = 82 
and it crosses the continuum at A^ = 100 for the nu- 
cleus to ^''"Zr. In this region the chemical potential is 
very small and almost parallel to the continuum limit. 
This means that the additional neutrons are added with 
a very small, nearly vanishing binding energy at the edge 
of the continuum. The total binding energies E for the 
isotopes above ^^^Zr are therefore almost identical. This 
has been recognized already in Ref. Q in RMF calcu- 
lations using the BCS approximation and an expansion 
in an oscillator basis, which is definitely not reliable for 
chemical potentials so close to the continuum limit. In 
the present investigation we obtain the same result treat- 
ing now the continuum and its coupling to the bound 
orbits by pairing correlations properly within full rela- 
tivistic HB-theory in coordinate space. 

To understand the kink in the neutron rms-radii shown 
in Fig. 1 we present in Fig. 2 the contributions {r)ij of 
the different quantum numbers to these quantities as a 
function of the mass number. It is clearly seen that the 
radii for the negative parity levels close to the continuum 
limit are responsible for the rapid increase of the neutron 
radius. Neutrons above the closed neutron core A^ = 82 
are filled into these orbits. As more and more neutron 
are added, 8^3/2 and 2^:7/2 ( after A^ > 88 ), 3pi/2 ( 
after A^ > 92 ) respectively become weakly bound, then 
the contribution of further continuum 2/5/2 and I/19/2 
become more and more important. Going from A = 122 
to A — 140 we observe an almost constant contribution 
of all the channels to the total rms matter radius except 
a sudden increase in the contribution of the 8^3/2, 2/7/2, 
8pi/2 a.nd 2^5/2 channels. This means that the giant halo 
in 124- 140 2j. formed by the occupation of all these 
levels in the respective nucleus. 
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In the upper panel of Fig. 3 we show the correspond- 
ing density distribution for neutrons and protons in the 
nucleus ^^^Zr. By dashed lines we show calculations for 
different values of the box size R =15, 20, 25, and 30 fm. 
It is clear that we need very large box sizes to describe 
the halo properly. For R =30 fm the neutron density is 
reliably reproduced only up to r = 25 fm, where it has 
decreased to 10~^ fm~^. The full line in the upper panel 
of Fig. 3 is an asymptotic extension to infinite box size. 
On the other hand, the density distribution inside the 
nucleus is reproduced properly even for small values of 
the box size. 

In the lower panel of Fig. 3 we show the relative con- 
tributions pnij of the different orbits characterized by the 
quantum numbers nlj with respect to the total neutron 
density pn- For comparison we also show the total neu- 
tron density in the shaded area in arbitrary units. As we 
see the halo is formed essentially by contributions from 
three orbits 3p3/2, 3pi/2, and 2/7/2- The most inner part 
of the halo (7 < r < 9 fm) the 2/7/2 orbit plays the dom- 
inant role. As can be seen in the lower part of Fig. 1, 
this orbit is slightly below the chemical potential and 
to the continuum limit in this nucleus. Further outside 
(10 < r < 15 fm) its relative contribution is strongly 
reduced because of the larger centrifugal barrier felt by 
the I = 3 orbit. In this region the orbit 3p3/2, which 
has nearly the same position as the 2/7/2 orbit, takes 
over. Because of the smaller orbital angular momentum, 
this orbit feels a reduced centrifugal barrier. For even 
larger distances from the center (r > 15 fm) its relative 
contribution is somewhat reduced and the 3pi/2 orbit 
gains importance. The 3^3/2 and the levels feel 

the same centrifugal barrier, but the latter is situated 
directly at the continuum limit and therefore it is more 
loosely bound than the other two orbits. 

In Fig. 4 we show for all the Zr-isotopes between 
A = 108 and A = 140 the occupation probabilities in the 
canonical basis of all the neutron levels near the Fermi 
surface, i.e. in the interval -20 < < 10 MeV. The 
chemical potential is indicated by a vertical line. For the 
mass numbers A < 122 the chemical potential lies sev- 
eral MeV below the continuum limit {E = 0) and there 
is only very little occupation in the continuum {E > 0). 
The nucleus ^^^Zr has a magic neutron number and no 
pairing. As the neutron number goes beyond this closed 
core, the occupation of the continuum becomes more and 
more important. Adding up the occupation probabilities 
v"^ for the levels with E > we find in the continuum 2 
particles for = 84, 4 for N = 86, 6 N = 88, roughly 3 
for N = 90, roughly 4 for A = 92, roughly 3 for A = 94, 
roughly 4 for AT = 96, roughly 5 for A" = 98, and roughly 
6 for A' = 100 (where the neutron drip line is reached). 

Summarizing our investigations, we predict neutron 
halos in the Zr nuclei close to the neutron drip line, which 
are composed not only by one or two neutrons, as is the 
case in the halos investigated so far in light p-shcll nuclei, 
but which contain up to 6 neutrons. This is a new phe- 
nomenon, which has not been observed experimentally 



so far. It would allow the study of collective phenom- 
ena in neutron matter of low density. This prediction 
is based on relativistic Hartree-Bogoliubov theory in the 
continuum. It combines the advantages of a proper de- 
scription of the spin orbit term with those of full Hartree- 
Bogoliubov theory in the continuum, which allows in the 
canonical basis the scattering of Cooper pairs to low ly- 
ing resonances in the continuum. A density dependent 
force of zero range has been used in the pairing chan- 
nel. It contains no free parameter, because its strength 
is adjusted for the isotope ^^^Zr to a similar calculation 
with Gogny's force DIS in the pairing channel. The halos 
are formed by two to six neutrons scattered as Cooper- 
pairs mainly to the levels 3p3/2, 2/7/2, ''^Pi/2-, and 2^^5/2- 
This is made possible by the fact that these resonances in 
the continuum come down very close to the Fermi level in 
these nuclei and by their coupling with the loosely bound 
levels just below the continuum limit. 
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Figure Captions 

Fig. 1 Upper part: rms-radii for neutrons and protons 
in Zr isotopes close to the neutron drip line as a 
function of the mass number A. Lower part: single 
particle energies for neutrons in the canonical basis 
as a function of the mass number. The dashed line 
indicates the chemical potential. 

Fig. 2 Contributions for the different channels 

with the quantum numbers I and parity as a func- 
tion of the mass. The left part show the orbits with 
positive parity and the right part shows the levels 
with negative parity. 
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Fig. 3 Upper part: neutron and proton density distribu- 
tion in ^^^Zr. Dashed lines indicate calculations for 
different values of the box size R. Lower part: rela- 
tive contributions of the different orbits to the full 
neutron density as a function of the radius. The 
shaded area indicates the total neutron density in 
arbitrary units. 

Fig. 4 The occupation probabilities in the canonical ba- 
sis for various Zr isotopes with mass number A as 
a function of the single particle energy. The chem- 
ical potential is indicated by a vertical line. For 
A > 122 we also show the number N^, of neutrons 
in the halo. 
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